CMSSM Spectroscopy in light of PAMELA and ATIC 
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Dark matter neutralinos in the constrained minimal supersymmetric model (CMSSM) may account for the 
recent cosmic ray electron and positron observations reported by the PAMELA and ATIC experiments either 
through self annihilation or via decay. However, to achieve this, both scenarios require new physics beyond 
the 'standard' CMSSM, and a unified explanation of the two experiments suggests a neutralino mass of order 
700 GeV - 2 TeV. A relatively light neutralino with mass around 100 GeV (300 GeV) can accomodate the 
PAMELA but not the ATIC observations based on a model of annihilating (decaying) neutralinos. We study the 
implications of these scenarios for Higgs and sparticle spectroscopy in the CMSSM and highlight some bench- 
mark points. An estimate of neutrino flux expected from the annihilating and decaying neutralino scenarios is 
provided. 

PACS numbers: 95.35.+d, 98.62.Gq, 98.70.Vc, 95.85.Ry 



I. INTRODUCTION 

It is generally accepted that nearly 23% of the universe's en- 
ergy density resides in the form of non-luminous 'dark' mat- 
ter iQl]. This is a new form of matter which is non-baryonic 
and manifests itself primarily through its gravitational interac- 
tions. The highly successful Standard Model (SM) of strong, 
weak and electromagnetic interactions does not possess a vi- 
able dark matter candidate. Thus, new physics beyond the SM 
is required to incorporate dark matter, and many potential dark 
matter candidates have been proposed in the literature \2\. 

Supersymmetry, more precisely MSSM, (minimal super- 
symmetric SM), with R-parity conservation, is arguably the 
most compelling extension of the SM. The MSSM predicts the 
existence of a stable new elementary particle called the neu- 
tralino (lightest supersymmetric particle). With mass of order 
100 GeV - TeV, the thermal relic abundance of the lightest 
neutralino has the right order of magnitude to account for the 
observed dark matter density. 

Many recent investigations of the MSSM have focused on 
a theoretically well motivated special case called the CMSSM 
dll] (constrained MSSM, based on supergravity) which is far 
more predictive than the generic MSSM version. The latter 
can have more than a hundred free parameters, in contrast to 
the CMSSM with just 5 or so parameters. In these investi- 
gations the stable neutralino is usually found to be relatively 
light (100-few hundred GeV), and its indirect discovery relies 
on detecting cosmic ray signals (including positrons, antipro- 
tons, gamma rays, etc.) from neutralino decay or pair anni- 
hilation in the galactic halo, galactic center, and the haloes of 
nearby galaxies. 

The PAMELA experiment is currently taking data of high 
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energy anti-proton and positron fluxes, and their most recent 
publication claims a significant positron 'excess' |3l with no 
corresponding anti-proton excess This result appears to 
confirm previous results from HEAT [6] and AMS |7] within 
the error bars. It has been pointed out that pulsars and/or other 
nearby astrophysical sources may account for the PAMELA 
results 1 8, 9]. More recently, the ATIC experiment |10] (see 
also PPB-BETS [11]) has reported an appreciable flux of elec- 
trons and positrons at energies around 100 - 800 GeV, which 
appears to be considerably higher than the expected back- 
ground at these energies. 

A unified explanation of the PAMELA and ATIC experi- 
ments involving neutralino as the dark matter candidate could 
be based on one of the following two mechanisms: 

• The dark matter is a stable neutralino with mass 
around 700 GeV which primarily annihilates into lep- 
tons through new interactions which lie outside the 
MSSM framework fl^j. Depending on the framework 
chosen, this scenario also invokes some 'boost' factor 
physics such as Sommerfeld enhancement ifisll . 

• The dark matter is not entirely stable ifTil [Tsll but ex- 
tremely long-lived, with a lifetime ~ 10^^ sec. For neu- 
tralino dark matter, one could introduce suitably 'tiny' 
R (or 'matter') - parity violating couplings which sat- 
isfy the lifetime constraint and allow the neutralino to 
decay primarily into leptons. The tiny (^ lO^^'') R- 
parity violating couplings can be understood through 
non-renormalizable couplings with additional discrete 
symmetries fl^. A simple example of this is pro- 
vided by the R-parity violating superpotential coupling 
LLE'^, which leads to a three-body decay mode for the 
neutralino ifTil] . With a neutralino mass of around 2 
TeV, this can simultaneously explain the PAMELA and 
ATIC data. 

If one ignores say the ATIC result, it is possible to explain 
the PAMELA observations with a decaying neutralino of mass 
~ 300 GeV, which would make supersymmetry, and in par- 
ticular the sparticles, far more accessible at the LHC. Even 
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lighter ('^ 100 GeV) neutralino mass is feasible in the annihi- 
lation scenario. 

Motivated by the PAMELA and ATIC observations we have 
performed an ISAJET |17] based analysis of CMSSM spec- 
trocopy in which particular attention is paid to those regions 
of the CMSSM parameter space which contain heavy (~ 300 
GeV - 2 TeV) neutralinos, with a relic abundance consistent 
with the WMAP 5 dark matter bounds 1 1], and which provide 
a unified explanation of the PAMELA and ATIC observations. 

The plan of the paper is as follows. We summarize the ob- 
servations by PAMELA/ ATIC experiments in Section II and 
review their implications in the context of dark matter de- 
cays or annihilations in Section III. In Section IV an estimate 
of neutrino flux expected from the annihilating and decaying 
neutralino scenarios is provided. Section V discusses how the 
CMSSM could be supplemented by new physics while pre- 
serving the neutralino cold dark matter framework. We outline 
in Section [Vl] the procedure that we use to scan the CMSSM 
parameter space and the various experimental bounds that we 
take into account. In Section IVIII we present plots display- 
ing the relevant CMSSM parameter space and highlight in Ta- 
bles U and a few benchmark points that are consistent with 
the PAMELA and ATIC measurements. Our conclusions are 
summarized in Section lVIIII 



II. OBSERVATIONS 

The spectra of cosmic ray electrons and positrons should 
have contributions from known sources including e~ acceler- 
ated in supernova remnants and from collisions between 
cosmic rays and interstellar protons [18]. Besides well known 
guaranteed backgrounds, any evidence for an additional com- 
ponent may carry indications of a new phenomon. Recently, 
the PAMELA |4] and ATIC |10] experiments have observed 
cosmic ray positrons and/or electrons within the energy range 
from several GeV up to a few TeV, and provided new hints for 
the long suspected excess in their fluxes f^J^ with respect to 
the model expectations |19]. The revitalized interest in this 
long standing puzzle has recently generated much interest in 
finding both astrophysical and perhaps more exotic explana- 
tions. 

The ATIC experiment has measured the combined electron 
and positron flux from ^30 GeV reaching energies up to sev- 
eral TeV as shown in Fig. 1. The primary electron spectrum 
is based on calculations of Ref. |19] (light shaded) is also 
shown. As the normalization of the primary electron spectrum 
is already quite uncertain, we renormalize this down by a fac- 
tor ~ 0.8 in order to match the low energy part of the ATIC 
data. Unlike electrons, while there are no primary sources of 
positrons in this model (see also (idH ). the expected secondary 
positron spectrum mainly due to cosmic ray interactions (dark 
shaded) is also presented. There is also a much smaller sec- 
ondary electron component (not shown in the figure). 

The PAMELA data is originally reported as the positron 
fraction of total electron/positron flux up to ^ 100 GeV |4]. 
In order to show both electron and positron data on the same 
figure, we have converted the reported ratios into absolute 



flux units using the primary electron spectrum. As shown 
in the figure, the corresponding positron flux based on the 
PAMELA data suggests a strong deviation from the secondary 
positron background. While the measurements by PAMELA 
below ^10 GeV might depend on solar modulation effects, 
the excess in the positron flux measurements with respect to 
the model expectations is evident at higher energies. 

The ATIC measurement of the total electron/positron flux 
also shows some excess beyond ^100 GeV, most prominently 
at ~ 700 GeV (once the baseline for the primary electron 
spectrum is normalized to the ATIC data at low energies where 
the positron contribution is small). Even though the reported 
excess in the ATIC data is not large compared to their uncer- 
tainties ijTod and possible systematic effects, and earlier ex- 
periments in a similar energy range show a rather large scatter 
(see e.g. compilation in [11]), ATIC observations nonetheless 
provide avenues for exploring new physics, especially when 
taken together with PAMELA. 

One can interpret the combined observations of PAMELA 
and ATIC as a hint of an unknown source capable of produc- 
ing comparable numbers of electrons and positrons in the en- 
ergy range of 10-700 GeV. For instance, if the flux of excess 
pairs follows, approximately, a power law of the form 

J{E) ~ 0.4 GeV^i sv'^ , (1) 

with a cut-off around ~ 700 GeV (solid line in Fig. 1), then 
both ATIC (dashed Une) and PAMELA (dot-dashed line) re- 
quirements can be reasonably satisfied. 

A very promising possibility that can account for the ex- 
cess is astrophysical in origin, i.e., pulsar wind nebulae which 
could accelerate pairs produced around a neutron star, 
and which would naturally yield equal numbers of positrons 
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FIG. 1: Cosmic ray electron (light shaded) and positron (dark 
shaded) spectra from the model of Ref. 1,19,1 . Data has been taken 
from Ref. IToll (squares), and Ref. 0] (circles). An additional 
contribution of pairs with spectra J{E) ~ 0.4 (_E/GeV)"^'^ 
GeV~^m~'^s~^sr~^ with a cut-off around ~ 700 GeV (solid line) 
can provide a better fit to both ATIC (dashed lined) and PAMELA 
(dot dashed line) measurements. 
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FIG. 2: Models of cosmic ray electron/positron background spectra and observations as in Fig.[T] Here the additional contribution (solid 
line) is assumed to be originated from dark matter decay (top panels) or annihilation (bottom panels) and aimed at explaining either only 
PAMELA (left panels; as a lower bound on dark matter mass) or also ATIC (right panels; as an upper bound on dark matter mass). See text 
for details. 



and electrons. Indeed, the tentative detection of TeV gamma 
rays from Geminga [21], a nearby pulsar, immediately sug- 
gests the possibility of particle accelaration up to multi-TeV 
energies (st) (provided that the origin of the gamma rays is 
inverse Compton scattering). Unlike electrons, which can 
also be accelarated in supernova remnants, positrons have no 
other known primary sources, and thus the resulting excess 
in positron flux could be very prominent and similar to the 
PAMELA observations I 



III. POSITRONS/ELECTRONS FROM DARK MATTER 

Another exciting explanation is that some exotic mecha- 
nism such as dark matter decay or annihilation, which are the 
main subjects of this study, might be responsible for the ex- 
cesses seen in observations. We will assume throughout that 
the dark matter is the lightest CMSSM neutralino. For a study 
on the signals from dark matter annihilation/decay products, 
an understanding of the distribution of dark matter within the 
galaxy is essential. Simulations of cold dark matter suggest 
that dark matter assembles in halos. The dark matter distribu- 
tion within the halo is usually taken to be spherically symmet- 



ric and the density profile can be expressed as function of the 
distance from the center of the halo. 



p{r) 



Po 



(r/rs)T[l + (r/rs)«]('5-T)/"^ 



(2) 



where the parameters (a, (3,j,rs) ~ (1, 3, 1, 20 kpc) for the 
NFW profile l'22'j, which also suggests that the density of dark 
matter at the solar circle is p{Rsc — 8.5 kpc ) = psc — 0.3 
GeV/cm^. 

As our main motivation is to account for the observations 
of ATIC and PAMELA, we shall focus on generic leptophilic 
scenarios in which the dominant annihilation/decay products 
are assumed to be positrons and electrons. In such models, 
once electrons/positrons are created, their density n at a given 
place, time and energy, is governed by the diffusion equation. 



— = V • [V{E)Wn] + — [£{E) n] + Q{r, E) 



(3) 



where Q{y, E) is the source term describing the particle in- 
jection rate as discussed below, 1{E) = ^o(-E'/GeV)~^ is the 
energy loss rate and 'D{E) = V^iE /GcYY is the diffusion 
coefficient, both assumed to be independent of space. 

The spectrum of a particle per dark matter annihilation, 
XX ~^ e+e", can be written as 4)°'{E) = 5{E — m^) in 
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terms of the Dirac 5 function, where is the mass of the 
dark matter particle. The corresponding source term is then 



>(r)' 
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2 TO^ 



r{E), (4) 



where the factor of 1 /2 arises by assuming that the dark mat- 
ter candidate is its own antiparticle. Here the term fB{(Jv) 
is the product of the annihilation cross section and the rela- 
tive velocity together with the requried 'boost factor fe' that 
may come from either new particle physics or astrophysical 
enhancements. 

If the dark matter decays into three light particles such as 
X e^e^v, the spectrum of each daughter particle, to a 
rough approximation, is given by cj)'^{E) ~ S{E — rn^/S) 
with the corresponding source term 



Q'(r,i?) = 



Psc 



Ps 



(5) 



where r is the lifetime corresponding to this decay channel, 
which may arise from new physics. Note that the prefactor 
Psc is introduced to make the term in the second paranthesis 
dimensionless and will cancel out in our final results. 

With these definitions of the source terms, we follow an ap- 
proach similar to Ref. [23, 24] in evaluating the local particle 
flux. Assuming steady-state conditions, typical in assessing 
the particle flux from dark matter annihilation/decay, the so- 
lution of the diffusion equation can be cast as 



n{E) 



dE (f){E )C{E,E 



(6) 



The coefficient k encodes particle physics input such that for 
annihilations k°' = i(pg^/m^)fB((Tw), while for decays k'^ ~ 
Psc/ (w^r). The spectrum of each given daughter particle per 
annhilation or decay is 0. The function C encodes dependence 
of the solution on the chosen dark matter halo profile as well 
as particle propagation and diffusion and has a form similar 
to / in Ref. [23] for annihilations up to multiplicative factors. 
The flux of particles is simply J{E) = (c/47r)n(£'), where c 
is the speed of light. 

The meaning of Eqn. |6]is as follows. For a particle to be 
observed locally at a given energy E, it has to be produced at 
a higher energy i? at a more distant location in the galaxy. 
Then the particle will lose energy and diffuse until it reaches 
the solar system, abiding by the constraints of the diffusion 
equation. As high energy particles loose energy fast, their ob- 
served flux should be dominated by nearby processes within 
a few kiloparsecs. Higher production rate for rather cuspy 
halo profiles, especially for annihilations around the central 
parts of the galaxy are only important at lower energies. Since 
we focus on leptophilic scenarios, plausible severe constraints 
due to final states other than pairs, such as antiprotons, are 
akeady avoided. Thus the exact parametrization of the under- 
lying dark matter halo profile is less relevant as long as the lo- 
cal normalization of dark matter density, psc — 0.3 GeV/cm'^ 
is consistent. 

Besides electrons/positrons, if other charged leptonic inter- 
mediate states exist, such as tau that decay to e^, the fit of the 



theoretical expectations to the observations might be visually 
improved. However, an exact reproduction of the observations 
cannot be essential since (1) the observations are highly un- 
certain, and more importantly (2) the background fluxes are 
not well known and already strongly dependent on models. 
Since an attempt to over-constrain the theories while assum- 
ing a perfect knowledge of the background electron/positron 
fluxes may not yield additional useful results, we choose to 
probe only the most generic features of the observations by 
R\MELA and ATIC and present our results for a NFW dark 
matter profile together with the MED diffusion model ||23|] 
with appropriate parameters {Dq = 0.0112 pc^/yr, (5 = 
0.7,^0 = 10"^^ GeV/s). 

Next we outline how the R\MELA and ATIC observations 
could be accounted for by a neutralino dark matter scenario 
through annihilation or decay, and what would be the ex- 
pected ranges of parameters such as dark matter mass, life- 
time or annihilation cross section under these assumptions. 
We identify two scenarios which can reasonably account for 
the positron excess observed by PAMELA with a relatively 
light (~ 100 - 300 GeV) dark matter candidate. In addi- 
tion we consider two scenarios with an appropriately heavy 
700 GeV — 2 TeV) dark matter candidate which provide a 
unified explanation of both the PAMELA and ATIC observa- 
tions. 

In Fig. 2, the additional contribution (solid line) is as- 
sumed to originate from dark matter decay or annihilation. 
The dashed line shows the expectation for total electron and 
positron flux when both e~ and e+ fluxes are added to the 
primary electron background (light shaded). This can be 
compared to the ATIC data which corresponds to the to- 
tal lepton spectrum. We also show the total positron spec- 
trum (dot-dashed line) by adding only the additional positron 
flux from dark matter annihilation or decay to the secondary 
positron spectrum (dark shaded), which is then compared to 
the PAMELA data. The theoretical fits are slightly smoothed 
out to simulate finite resolution of experiments. The four pan- 
els displayed are: 

• Top-Left: A decaying dark matter scenario with a life- 
time r ~ 2 X 10^^ s and ^ 0.3 TeV which is the 
minimal mass required to account for PAMELA (dot- 
dashed line) without causing significant tension with 
ATIC (dashed-Une). 

• Top-Right: A decaying dark matter scenario with a life- 
time T ~ 3 X 10^^ s and ~ 2 TeV which is the 
maximum mass that can account for ATIC observation 
with rather sharp cut-off around ^ 0.7 TeV (dashed 
line) and also account for PAMELA (dot-dashed line). 
Thus, masses in the range of ^ 0.3-2 TeV with a life- 
time of ^ 10^^ s can be regarded as acceptable for a 
decaying neutralino dark matter scenario. 

• Bottom-Left: An annihilating dark matter scenario with 
fB(cr?;) - 1.5 X 10^25 ^j^3/g ^jjjj _ 0.1 TeV, 
which corresponds to the minimal mass required to ex- 
plain PAMELA (dot-dashed line) without significantly 
exceeding the ATIC measurement (dashed-line). 
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• Bottom-Right: An annihilating dark matter scenario 
with fsiav) - 6 X 10^24 ^^3f^ ^ q 7 j^y 

that can simultaneously fit the PAMELA excess (dot- 
dashed line) and ATIC excess and also the observation 
of rather sharp cut-off around ^ 0.7 TeV (dashed line). 
For an annihilating scenario, the preferred mass ranges 
are - 0.1-0.7 TeV with iB{crv) - 10^24 cm^/s. 

The predictions for dark matter annihilation and decay 
barely differ only at low energies as can be seen in Fig. 2, 
where the flux from the former is larger This is due to the fact 
that the annihilation rate is proportional to the square of the 
dark matter density and any enhancement towards the Galac- 
tic Center will cause an enhanced contributions at lower ener- 
gies. 



IV. NEUTRINOS FROM DARK MATTER 

An interesting application of these results is the plausible 
presence of other leptonic final states such as neutrinos. Dark 
matter annihilating dominantly into neutrinos has been con- 
sidered in Ref. luSl l26ll and more recently in the context of 
more recent observations 28]. For instance, neutralino 
decay may yield neutrinos besides pairs of e^. Similarly, it 
is conceivable that in annihilation, neutrinos and pairs are 
produced in comparable numbers. 

To further illustrate the implication of the four scenarios, 
we calculate the total neutrino flux (utilizing prescriptions 
outlined in Refs. |2^ |2^) using the above parameters and 
present the results in Fig. 3. The total neutrino flux for these 
bumps is 6 X 10~^ cm~2s~^sr~^ when averaged over the 
whole sky (both smoothed with a 10% Gaussian in the figure). 
The fluxes from the decay (thick lines) and annihilation (thin 
lines) scenarios are comparable as their main normahzations 




1 2 3 4 

10 10 10 10 

E [ GeV ] 

FIG. 3: The atmospheric neutrino spectrum averaged over whole sky 
(solid line) is compared to the total neutrino flux associated with the 
parameters used in top-left ( top-right, bottom-left, bottom-right) 
panel of Fig. 2 which is denoted with a thin-dot-dashed (thin-dashed, 
thick-dot-dashed, thick-dashed) line. See text for details. 



are provided by the requirement to account for the R^MELA 
and ATIC observations. The fluxes are also averaged over the 
whole sky which dampens any relative enhancement of anni- 
hilation compared to decay scenarios, especially towards the 
center of the Galaxy. 

These can be compared to the atmospheric neutrino 
spectrum (solid line) which is based on measurements by 
Refs. |30, 31, 32, 33^^34] and also agrees with the theoretical 
modelling of Ref. i35ll . The bumps at ~ 100 GeV which are 
devised to account for PAMELA only fall short by several or- 
ders of magnitude below the atmospheric neutrino flux. How- 
ever at energies ^ 700 GeV, the coiTesponding neutrino fluxes 
are within reach of the atmospheric neutrino background, thus 
offer hope to test this scenario in the near future. 

As we have focused here on leptophilic scenarios, we do not 
discuss models in which gamma ray fluxes are among the pri- 
mary decay or annihilation products (see e.g. Refs. ll36l IstIi 
where the dominant dark matter annihilation/decay products 
are gamma rays). However, gamma rays might accom- 
pany the production of pairs such as through internal 
bremsstrahlung |38]. Also inverse Compton and synchrotron 
emission might be produced from the electrons and positrons. 
In such cases, especially for annihilations, dark matter profiles 
that are shallower around the center of the Galaxy compared 
to the NEW parametrization might be preferred Il39tl . 



V. FROM OBSERVATIONS TO CMSSM 

In order to understand the unexpectedly high electronic flux 
observed by PAMELA and ATIC the CMSSM must be sup- 
plemented by new physics. If one wishes to preserve the basic 
neutralino cold dark matter framework, which we propose to 
do in this paper, at least two options are available. 

Perhaps the simplest modification one could contemplate is 
to include R-parity violating superpotential coupling(s) with 
suitably tiny dimensionless coefficient(s), such that the neu- 
tralino primarily decays into leptons, with lifetime ^ lO^^sec. 
A simple example is provided by the coupling LLE'^ fTil], 
which leads to 3-body leptonic decays of the neutralino. This 
modification has an important advantage in that it preserves 
the basic structure of the CMSSM, especially by leaving intact 
the standard calculations for estimating the relic neutralino 
abundance. The implications for CMSSM spectroscopy for 
this case which has a rather heavy neutralino (^ 2 TcV) are 
discussed in Section lVTIl 

Dark matter neutralinos annihilating in the halos of galax- 
ies provide an exciting new source for cosmic rays including 
positrons. To explain the PAMELA and ATIC observations in 
terms of neutralino annihilations, the CMSSM must be sup- 
plemented by new physics which should accomplish the fol- 
lowing three things: allow neutralino annihilations primarily 
into leptonic channels, produce a boost factor of order 10'' or 
so through Sommerfeld enhancement lfl3ll . and last but by no 
means least, ensure that the relic abundance of neutralinos is 
consistent with WMAP bounds. Several interesting propos- 
als for annihilating WIMPS that explain PAMELA and ATIC 
have been put forward 1 12]. However, it is not obvious if they 
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FIG. 4: Plots in (mi/2, mo) and (mi/2, tan/3) planes for /i > (left panel) and fi < (right panel). Gray points satisfy constraints from 
colliders ( BR{Bs — > jj,'^ fj,~), BR{B ~~* Xs')), and the chargino and Higgs mass bounds). Light blue points satisfy the WMAP upper bound 
on dark matter relic abundance. Green {m.o < 0.3 TeV), red (0.65 TeV < m^o < 0.75 TeV), black (1.9 TeV < m^o < 2.1 TeV), and 

Ai Xi Xi 

dark blue ("i^o < 2.5 TeV) points satisfy both the upper and lower bounds on dark matter relic abundance. Green, red and black points are 
subsets of dark blue points. Dense regions in all figures correspond to parameter values that are especially relevant for PAMELA and ATIC, 
and for which additional data has been accumulated. 



can be easily adapted to the CMSSM scenario we are consid- 
ering in which the neutralino makes up the observed cold dark 
matter in the universe. Thus, even though an explicit particle 
physics construction for neutralino annihilation satisfying the 
above conditions is not yet available, we will in the following 
discussion present spectroscopy results which cover a wide 
range of neutralino masses currently favored by PAMELA and 
ATIC, including the 700 GeV or so value motivated from neu- 
tralino annihilation. 



VI. PHENOMENOLOGICAL CONSTRAINTS AND 
SCANNING PROCEDURE 

We employ ISAJET 7.78 package ITtIi to perform random 
scans over the parameter space. In this package, the weak 
scale values of gauge and third generation Yukawa couplings 
are evolved to A/qut via the MSSM renormalization group 
equations (RGEs) in the DR regularization scheme, where 
Mqut is defined to be the scale at which gi = g2- We 
do not enforce an exact unification of the strong coupling 
93 = 91 — 92 at Mqu'y, since a few percent deviation from 
unification can be assigned to unknown GUT-scale threshold 
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FIG. 5: Plots in (mo, ^o) and (mi/2, Aq) planes for > (left panel) and jj, < (right panel). Color coding same as in Fig.U 



corrections f4^. At A/qut, the CMSSM boundary condi- 
tions are imposed and all the SSB parameters, along with 
the gauge and Yukawa couplings, are evolved back to the 
weak scale Mz- In the evaluation of Yukawa couplings the 
SUSY threshold corrections |41] are taken into account at 
the common scale A/susy = ^/^h ■ The entire pa- 
rameter set is iteratively run between Mz and Mqut using 
the full 2-loop RGBs until a stable solution is obtained. To 
better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and 
the SSB parameters rrii are extracted from RGBs at multiple 
scales nii — m.i{mi). The RGB-improved 1-loop effective 
potential is minimized at an optimized scale A/susy, which 
effectively accounts for the leading 2-loop corrections. Bull 
1-loop radiative corrections are incorporated for all sparticle 
masses. 



The requirement of radiative electroweak symmetry break- 
ing (RBWSB) |42] puts an important theoretical constraint 
on the parameter space. Another important constraint comes 
from limits on the cosmological abundance of stable charged 
particles |43l] . This excludes regions in the parameter space 
where charged SUSY particles, such as fi or ti, become the 
lightest supersymmetric particle (LSP). We accept only those 
solutions for which the neutralino is the LSR 

We have performed random scans for the following param- 
eter range: 

< mo < 20 TeV, 
0< TOi/2 < 10 TeV, 
-3< Ao/mo <3, 
5< tan/3 < 58, (7) 
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FIG. 6: Plots in (rriA, rn^o) and (m-±, m-o) planes for u > (left panel) and u < (right panel). Gray points satisfy constraints from 

colliders ( BR{Bs fi'^ fi^), BR{B -^^^7), and the chargino and Higgs mass bounds). Light blue points satisfy the WMAP upper bound 
on dark matter relic abundance. Dark blue points satisfy both the upper and lower bounds on dark matter relic abundance. Orange and brown 
points satisfy the constraint from Aa^ (for /i > 0, as expected). Orange points satisfy only the lower bound on dark matter relic density, 
while brown ones satisfy both the upper and lower bounds. Vertical lines correspond to m^o = 0.3, 0.65, 0.75, 1.9 and 2.1 TeV. ATIC 
imposes 0.65 TeV < m^o < 0.75 TeV (1.9 TeV < m^o < 2.1 TeV) based on annihilating (decaying) neutralinos. PAMELA requires 
— TeV (m^o > 0.3 TeV) based on annihilating (decaying) neutralinos. 



with ^ > and < 0, and rrit = 172.6 GeV Q. 

After collecting the data, we use the IsaTools package ||45|] 

to implement the following phenomenological constraints: 
m-± (chargino mass) > 103.5 GeV I 43ll . 
nih (lightest Higgs mass) > 114.4 GeV jiell . 
BR{Bs fi+n') < 5.8 X 10-* S, 
2.85 X 10-4 < BR{B Xs"/) < 4.24 x lO"'^ {2a) B, 
^CBuh^ = 0.11ltS;:;jl^ {2a) H, 
3.4 x 10^1° < Aa^ < 55.6 x IQ-io {3a) IH]. 
We have applied the constraints from experimental data 



successively on the data that we acquired from ISAJET. As 
a first step we apply the constraints from BR{Bs — > fi^n^), 
BR{B -Vs7), chargino mass, and Higgs mass. We then 
apply the WMAP upper bound on the relic density of cold 
dark matter followed by the constraint on the muon anoma- 
lous magnetic moment = {g — 2)^/2 at the 2>a allowed 
region. Finally, we apply the lower bound on the dark matter 
relic abundance. The data is then plotted showing the succes- 
sive application of each of these constraints. 



9 



jf/ > 



fi<0 






3 




2.5 




2 












1.5 


1 






1 




0.5 









0.5 1 1.5 

m^o{TeV) 



2.5 



10 
8 
6 
4 
2 




I- 



0.5 



3 

2.5 P 

2 : 
1.5 

1 

0.5 




1 1.5 

m^o{TeV) 

fi<0 
— I : 



Q QJ 1 1.5 

m^o [TeV] 



2.5 




2.5 



FIG. 7: Plots in (rrig, m^o ) and (mi, m^o ) planes for /i > (left panel) and fi < (right panel). Color coding same as in Fig. |6] Vertical lines 
correspond to m^o = 0.3, 0.65, 0.75, 1.9 and 2.1 TeV. ATIC imposes 0.65 TeV < m^o < 0.75 TeV (1.9 TeV < m^o < 2.1 TeV) based 
on annihilating (decaying) neutralinos. PAMELA requires m^o > 0.1 TeV (m^o > 0.3 TeV) based on annihilating (decaying) neutralinos. 



VII. IMPLICATIONS FOR CMSSM 

PAMELA and ATIC, as discussed previously, have set con- 
straints on the mass and other properties of the dark matter 
candidate which, in our CMSSM case, is the neutralino. The 
preferred mass values are as follows: 

• m^o > 0.3 TeV: This region is interesting if we wish 
to explain the PAMELA results based on decaying neu- 
tralinos. 

• 1.9 TeV < m^o < 2.1 TeV: This relatively naiTow 
range of neutralino mass explains the results from both 
PAMELA and ATIC based on decaying neutralinos. 



• m^o ^ 0.1 TeV: This region of neutralino masses 
can explain the PAMELA observation if we assume that 
neutralinos annihilate, rather than decay, preferentially 
into electrons and positrons. 

• 0.65 TeV m^o < 0.75 TeV: This region can ex- 
plain the results from both PAMELA and ATIC based 
on annihilating neutralinos. 

We now present the results of the random scan focusing at- 
tention on the above mentioned neutralino mass ranges. In 
Fig.|4]we plot the results in the (mi/2i "^o) and {mi/2, tan /?) 
planes for /i > (left panel) and /i < (right panel). All 
of these points satisfy the theoretical requirement of REWSB 
and correspond to a neutralino with mass less than 2.5 TeV. 
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FIG. 8: Plots in (mh, m^o) and (m?, m-o) planes for /I > (left panel) and ^ < (right panel). Color coding same as inFig.|6] Verticallines 
correspond to m^o = 0.3, 0.65, 0.75, 1.9 and 2.1 TeV. ATIC imposes 0.65 TeV < m^o < 0.75 TeV (1.9 TeV < m^o < 2.1 TeV) based 
on the model of annihilating (decaying) neutralinos. PAMELA requires m^o > 0.1 TeV (ni^o > 0.3 TeV) based on annihilating (decaying) 
neutralinos. 



In addition, these points satisfy the various experimental con- 
straints listed earlier Gray points satisfy the constraints from 
BR{B, fJ.^^J.^), BR{B Xs-i), and the Higgs and 
chargino masses. After application of the WMAP 5 upper 
bound on neutralino dark matter (light blue points), the al- 
lowed region is drastically shrunk, with the remaining gray 
points associated with an unacceptably high dark matter relic 
density. The other colors (green, red, dark blue and black) in 
Fig- HI correspond to different ranges of neutralino mass such 
that the neutralino satisfies both the WMAP 5 upper and lower 
bounds on dark matter relic abundance. The specific mass 
ranges have been picked out as interesting in trying to explain 
the data from PAMELA and ATIC. These are; 



• Green points correspond to m^o < 0.3 TeV. 

• Black points correspond to 1.9 TeV < m^o < 
2.1 TeV. 

• Red points correspond to 0.65 TeV < m^a < 
0.75 TeV. 

Note that if we want to explain just the PAMELA data through 
neutralino annihilation with a suitable 'boost' factor, there is 
almost no constraint on the CMSSM parameter space as it 
requires in this case m^a 0.1 TcV, which is satisfied by 
almost all CMSSM allowed points. 

The {mi/2, mo) plane for /i > is distinctly different from 
that of /i < 0. With /i < 0, the ratio between toi/2 and m^a 
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TABLE I: Sparticle and Higgs masses (in units of GeV), with rrit = 172.6 GeV and /i > 0. We present three benchmark points. Also 
included are the spin-independent neutralino-nucleon interaction cross-sections. Point 1 can explain results from PAMELA based on either a 
decaying or an annihilating neutralino model. Point 2 (Point 3) can explain ATIC and PAMELA based on the model of neutralino annihilation 
(decay). 



is fixed (more precisely in a narrow range), and the neutralino 
is always bino-like. While this is true also for the case /i > 
for low values of mi/2, we find solutions in which ^nd 
m^o lose this correlation for higher values, and the neutralino 
develops a strong higgsino component. This is the reason for 
the difference in the (jni/2, itiq) planes between /i > and 
/i < 0. This difference is more apparent in Fig. |7] where we 
plot results in the (nig, m^o) plane. 

In Fig.|5]we present the results in the (mo, Aq) and (mi/2, 
Aq) planes for both /i > (left panel) and /i < (right panel). 
Color coding is the same as in Fig. H] Even though it is not 
explicitly shown in Fig. |5] the constraint from Aa^ favors a 
small magnitude (< 4 TeV) for Aq and /i > 0. 

Let us now take a closer look at the sparticle spectroscopy 
that ai-ises in the CMSSM. In Figs. |6l Ul and |8] we present 
the results in the (m^, m-j^o), (w-±, m^o), (rrig, m^o), (toj, 
m^o), (m/i, m^o) and (mf, m^o) planes for both u > 
(left panels) and /i < (right panels). In these three figures 
we use consistent color coding. In particularl, points in gray 
satisfy the constraints from colliders on BR{Bs /i+/x"), 
BR{B Xsj), and the Higgs and chargino masses. Shown 
in light and dark blue are points that further satisfy bounds on 



the density of dark matter set by WMAP. Light blue points 
only satisfy the upper bound set by WMAP while dark blue 
ones satisfy both upper and lower bounds. Finally, we also 
show in orange and brown colors the regions allowed by the 
bounds on Aa^. Orange points satisfy only the lower bound 
on dark matter relic density while brown points saturate the 
WMAP bound on the dark matter relic abundance. 

The vertical lines in Figs.|6]|2]and[8]correspond to m^o — 
0.3,0.65,0.75, 1.9 and 2.1 TcV. These lines isolate the' in- 
teresting ranges of neutralino mass from the point of view of 
PAMELA and ATIC as discussed above. 

Let us focus on trying to explain both PAMELA and ATIC 
using a decaying neutralino. This corresponds to the interval 
between the vertical lines m^a = 1.9 and 2.1 TeV in Figs.|6] 
|7] and [8] Of particular interest is the t (stop) coannihilation 
region in Fig.|2l In Fig. [8] we see that the stau is quite heavy 
for 1.9 TeV < m^o 5: 2.1 TeV, so that coannihilation oc- 
curs mainly with the stop. The chargino can also be relatively 
'light' in these regions. However, it is clear from the (m^±, 
m^o) (Fig.|6]l and (m^, m^a) (Fig.|2ll planes that the chargino 
and the stop are not simultaneously 'light' . Therefore, to si- 
multaneously explain both PAMELA and ATIC using a decay- 
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TABLE 11: Sparticle and Higgs masses (in units of GeV), with mt — 172.6 GeV and /i < 0. We present three benchmark points. Also 
included are the spin-independent neutralino-nucleon interaction cross-sections. Point 1 can explain results from PAMELA based on either a 
decaying or an annihilating neutralino model. Point 2 (Point 3) can explain ATIC and PAMELA based on the model of neutralino annihilation 
(decay). 



ing neutralino requires a very heavy sparticle spectrum in the 
CMSSM, and the only particle that can be observed relatively 
easily at the LHC is the lightest Higgs boson. 

If we focus our attention on explaining only the PAMELA 
results (m^o > 0.3 TeV), we obtain a sparticle spectrum that 
is much more accessible at the LHC. As shown in Fig. |6]in 
this case we can find a relatively light chargino (~ 500GeV) 
and CP-odd and charged Higgs with masses ~ 600GeV. 

While we do not currently have a nice theoretical model for 
neutralino annihilations within the framework of CMSSM, it 
is plausible that such a model may be constructed in the future. 
If we allow for this possibility, it makes sense to explain ATIC 
and PAMELA in terms of neutralino annihilations and study 
the coiTesponding CMSSM sparticle spectroscopy. The ad- 
vantage here is that we need a significantly lighter neutralino 
to explain the ATIC and PAMELA observations, and the cor- 
responding spectroscopy is much more exciting from the point 
of view of the LHC. 

In Table HI (Table Hlb we present a few representative points 
for /i > (// < 0). All of these points correspond to a bino- 
like neutralino and have been chosen such that: 

• Point 1 may be used to explain PAMELA based on ei- 



ther decaying or annihilating neutralinos. 

• Point 2 may be used to explain both ATIC and 
PAMELA based on annihilating neutralinos. 

• Point 3 may be used to explain both ATIC and 
PAMELA based on decaying neutralinos. 

We also calculate for each benchmark point the neutralino- 
nucleon cross section. For /i > 0, the cross sections coiTe- 
sponding to the benchmark points 1 and 2 are a few orders of 
magnitude smaller then the exclusion limits set by the cuiTent 
experiments for direct dark matter detection such as CDMS 
I50i1 and XENON 10 [5 1.1 . The cross sections in the case of 
IJ, < are smaller still. 



Vm. CONCLUSIONS 

We have explored the implications for Higgs and sparticle 
spectroscopy of the CMSSM model in the light of recent cos- 
mic ray observations reported by the PAMELA and ATIC ex- 
periments. Our investigation is based on the premise that the 
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lightest CMSSM neutralino comprises the dark matter in the 
universe, and that the underlying physics should provide a uni- 
fied explanation, either through neutralino self annihilation or 
via its decay, for the positron and electron-positron excess re- 
ported by these two experiments. We have identified some 
benchmark points in the CMSSM parameter space which are 
consistent with this unified explanation. The corresponding 
sparticle mass spectra turns out, as expected, to be relatively 
heavy. The outlook for new particle discovery at the LHC is 
considerably enhanced if the ATIC data is simply ignored. In 
one example based on the decaying neutralino scenario, we 
find a relatively light chargino (~ 500 GeV) and CP-odd and 



charged Higgs with masses ~ 600 GeV. The CMSSM anal- 
ysis reported here can be extended to other related scenarios 
with neutralino LSP such as those based on SU(5) 1521 . 
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